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Sliding-Mode Control of Differential- 
Difference Systems 


Sliding modes in ordinary differential equa- 
tions. 


x(t) = Ax(t ) + Bu(t). 

(1) 

\ / if s(x) > 0 

if.(.)<0 

(2) 

"Sliding modes” in continuous-time 
ence systems. 

differ- 

x(i + T ) = Ax(t) + Bu{t). 

(3) 

S(x(t + r)) = 0 => u(t). 

(4) 
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We consider two configurations: 


Configuration A 

x(t) = Anx(t) + A\ 2 z(t) 
z(t + t ) = A 2 ix(t) + A 22 z(t) + B 0 u(t). 

Configuration B 

z(t + t) = A n z (t) + A 12 x(t) 

x(t) = A 2 ix(t) + A 22 z(t) -f Bou(t). 

It is assumed that x e R ni , z e R n 2 and 
u e R m . An, A 12 . A 2 i, A 2 2, B 0 are con- 
stant matrices of appropriate dimensions. 



Sliding-Mode Control Design 


A 12 = B X C 2 (5) 

Quasicontrol: 

v — C 2 z ( 6 ) 

Configuration A 

1. Sliding mode in differential subsystem. 

v — v*(x) => Sq(x ) = 0 (7) 

2. Sliding mode in difference subsystem. 

u => S(x, z ) = v*(x ) — C 2 z = 0 ( 8 ) 
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Configuration B 

1. “Sliding mode" in difference subsystem. 

v = Dz(t) => S Q (z ) = 0 (9) 

2. Sliding mode in differential subsystem. 


u S — Dz(t) — C' 2 x(f ) — 0 


( 10 ) 



Sliding Mode Control of Nondispersive 
Flexible Structures 


Flexible rod in compression. 


d 2 Q(t,x ) _ d 2 Q(t,x ) 
dt ^ dx* 

(id 

8Q(t,0) 

at = tt(t) 
dQ(t,l) _ 0 

dt 


Laplace transform approach. 


p 2 Q(p,x) = Q"(Pix) 

(12) 

Q'(p,o) = -w(p) 

(13) 

0 / (p,i) = o, 

(14) 


where Q(p,x) = CQ(t>x), w(p) = Cu(t). 



The solution of the boundary value problem 


Q(p,x) 


e p(x-l) _|_ e — p(x-l) ^ 
e p _ e -p p 


If an output variable is 


2/(0 = 1 ) 


then 

y{p) = <?(p, 1) = e p3 e ' -pp ^( p )~ 

In the time domain: 

2/U + 1) - y(t - 1) = 2u(t) 


or 


w(p)-( 15) 


(16) 


(17) 


(18) 


2/(0 - p(* - 2) = 2u(i - 1). 


(19) 



The equation can be written in the form of 
the difference-differential system as 


Configuration A: 

y(t) = zi(t) 
zi(t+l) = ^2(0 + 2u(t) 
z 2^ +1) = *l(0 


Configuration B: 

2/lO + l) = 1 / 2(0 + 2v(t) 
2 / 2 U+ 1 ) = 3/l(0 
i)(0 = u(0» 

where y x (i) = y(t). 
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Control Design 
Configuration A 


Zi (t) = -A sgn(y(t)). 

The equality is valid if 

-s(0 = ^i(0 + A sgn(y(t )) = 0. 

To achieve this the control should be 
u(0 = -^22(0 - ~A sgn(y(t + 1)). 

y(t + 1) = y(t) + j zi(r)dT. 


( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 
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As a result control is 

u(t) = ~^y(t - 1) - ~\sgn(y(t) + 

y(t - 1) - y(t - 2) -f 2 / u(t)cIt). 

Jt-l 

With this control the system is stabilized in 
finite time. 


Control Design 
Configuration B 

5(f) = (1 - A)y 2 (0 + 2v(f), (24) 

where |A| < 1. If the control is 

u(t) = —\xsgn (2 v(t) + (1 - A )j/(t - 1)) - 

(l-A)y(f- 1) 

then 

a — —2 nsgn(s). (25) 


342 



Rod with attached masses 


d 2 Q(t,x ) d 2 Q(t,x ) 

3t 2 dx 2 

(26) 

at = ~ u(t) 

(27) 

a«?(t,i)_ a 2 c?(t, i) 

^ dt 2 

(28) 

Configuration A: 


*l(0 = ^2(0 

*2(0 = -*2(0 + ^1 (0 

Zi(t) = Z 2 (t — 1) + 2u(t — 1) 



*2 (0 = -2l(t - 1) + 2x 2 (t - 1), 
where y(t ) = Q(t, 1) = xi(t). 

u(t) = -~z 2 (t) - ~fzsgn(Xxi(t + 1) + 
X2 (t + !))■ 
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+ 


State Estimation and Prediction 

Extrapolator 

xi(t + 1) 1 

x 2 (t + l)J U; 

<t>(f + 1 — T)/i(r)dr, 

where 
Mr) = 


i: 


Xl(t) 

X 2 (t) 


0 


zi(r — 2) + 2 u(t — 1 ) — 2x 2 {t — 2 ) 


Observer 

xi(t) = x 2 (t) + Li(xi(t) - i/(0) 

x 2 (i) = -x 2 (f) + ^i(t) + L 2 (xi(t) - y(0) 

Zl(t) = z 2 (^ — 1) + 2ti(< — 1) + 

L3(xi(t) - 2/(0) 

£ 2 (0 — -2i(t - 1) + 2x 2 (t - 1) + 

L 4 (xi(t) - y(t)). 


The same approach can be used for systems 
of connected rods with many attached masses, 
multiple controls and also for the case of dis- 
tributed actuators 


9 2 Q(t,x) 

dt 2 


d 2 Q(t , x ) 

dx 2 


+ <p(x)u(t) 


(29) 


where <^(x) is quasipolynomial. In all cases 
the solution of the boundary value problem 
for Laplace transformed variables leads to 

Configuration A or Configuration B. 


Dispersive Structures 


Euler-Bernoulli beam. 


d 2 Q(t,x) _ d 4 Q(t,x) 
dt 2 ~ dx 4 

(30) 

Q(t,0) = 0 

(31) 

Q' x (t. 0) = 0 

(32) 

Q'L(t, 1) = 0 

(33) 

Qxxx(l, 1) = «(»)■ 

(34) 


Second order dispersive structure. 


d 2 Q(t, x) d 2 Q(t, x) dQ(t,x) 

- = a(x) — H*>(x) 57 -( 35 ) 


dV 


dx" 


dx 
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General fourth order equation. 


d 2 Q(t,x) d 4 Q(t,x ) 

= -a(x) h 


dv 


dx 4 

b( + c (x) dQ(t ' x) 


The boundary conditions 


dx“ 


dx 


Q(f,0) = 0 

( 36 ) 

9Q(t,0) _ 

dx 

( 37 ) 

t) 2 Q(t, 1) 

0 = Wi(0 

dx 2 

( 38 ) 

d\Kt, 1) 

= u 2 (0- 

dx 3 

( 39 ) 

Integral Transform 


P(t,0=[ T>(Z,x)Q(t,x)dx 

Jo 

(40) 

If V satisfies an adjoint homogeneous bound- 

ary value problem then P(f,z) 
equation 

satisfies 

d 2 P(t,o a 2 P ((,0 , ,,, 

dt 2 di 2 

(41) 

AO = -a(0)2?U,0). 

(42) 
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Euler-Bernoulli beam 


If V is a solution of the boundary value prob- 

lem: 



d 2 V(i,x) _ 


(43) 

1 

1 

CM 

uy* 

cte 4 

0({,0) = 0 


(44) 

!>'({, 0) = 0 


(45) 

©"({, 1) = 0 


(46) 

©"'({, 1) = 0. 


(47) 


then P(t ,£) satisfies an equation: 


d 2 P(t,Z) 

Ot 2 


" »t 2 " 




where 


(48) 


tf>(0 = -©(?,!)• 


(49) 
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Initial values: V(0,x ) and D£(0,x). 

If Z>£(0,x) = 0 then 

P*(t, 0) = 0. (50) 

The possibility to choose £>(0,x) is an addi- 
tional degree of freedom that can be used to 
assign the desired value of <p(0- 

Nonsingularity condition 

P = 0 =i> Q = 0. (51) 

Output 

y (t) = P(t,0) = f 1 V(0,x)Q(t 1 x)dx. (52) 
JO 
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boundary 

controls 

u(t) 




aw 


Conclusions 

Sliding mode control became very popular recently be- 
cause it makes the closed loop system highly insensi- 
tive to external disturbances and parameter variations. 
Sliding algorithms for flexible structures have been 
used previously, but these were based on finite-dimensional 
models. An extension of this approach for differential- 
difference systems is obtained. That makes it possible 
to apply sliding-mode control algorithms to the vari- 
ety of nondispersive flexible structures which can be 
described as differential-difference systems. 

The main idea of using this technique for dispersive 
structures is to reduce the order of the controlled part 
of the system by applying an integral transformation. 

We can say that transformation “absorbs" the dis- 
persive properties of the flexible structure as the con- 
trolled part becomes dispersive. 
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